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1 . INTRODUCTION 

Several  authors  [1]  - [5]  have  discussed  methods  of  deriving  the 
equation  of  motion  of  a generalized  elastic  foundation,  which  often  is 
treated  as  a thin  elastic  layer  of  thickness  h that  has  one  face  bonded 
to  the  surface  of  a rigid  substrate.  Such  models,  which  include  the 
effects  of  foundation  inertia  and  shear  deformation,  provide  more 
realistic  descriptions  of  the  dynamic  response  of  elastic  foundations 
than  does  the  classical  Winkler  model,  which  treats  the  normal  stress 
at  the  surface  of  the  foundation  as  linearly  proportional  to  its  trans- 
verse deflection  w,  i.e., 

P2  = -kw  (1) 

However,  upon  examining  the  simplest  foundation  models  proposed  in 
references  [1],  [2],  and  [4]  in  the  case  of  slender  beams  on  elastic 
foundations,  it  becomes  evident  that  these  theories  do  not  yield  the 
correct  lowest  frequencies  for  the  depth-shear  and  depth-stretch  modes 
of  free  vibration.  (The  depth-shear  and  depth-stretch  modes  considered 


Tv.  Z.  Vlasov  and  U.  N.  Leont'ev.  Beams,  Plates  and  Shells  on  Elastic 
Foundations.  NASA  TT  F-357  (1966) 

^N.  S.  V.  Kameswara  Rao,  Y.  C.  Das,  and  M.  Anandakrishnan.  Variational 
approach  to  beams  on  elastic  foundations.  Proc.  Amer.  Soc.  Civ.  Engrs., 
J.  Engng.  Mech.  Div.,  97,  pp.  271-294  (1971). 

^E.  H.  Dowell,  Dynamic  analysis  of  an  elastic  plate  on  a thin,  elastic 
foundation.  J.  Sound  Vi b . , 35^,  pp.  343-360  (1974). 

^N.  S.  V.  Kameswara  Rao,  Y.  C.  Das,  and  M.  Anandakrishnan.  Dynamic 
response  of  beams  on  generalized  elastic  foundations.  Int.  J.  Solids 
Structures,  1_1_,  pp.  255-273  (1975). 

^G.  L.  Anderson.  The  influence  of  a Wieghardt  type  elastic  foundation 
on  the  stability  of  some  beams  subjected  to  distributed  tangential 
forces.  J.  Soi in  I 'M h . (to  be  published). 
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here  are  quite  analogous  to  the  well-known  thickness-shear  and  thick- 
ness-stretch modes  that  are  known  to  exist  in  infinite  elastic  plates. 

Related  or  analogous  mathematical  problems  arise  in  the  design  of 
ultrasonic  devices  which  contain  crystal  bars  and  plates  that  are 
rigidly  clamped  on  one  face  due  to  the  mode  of  mounting. 

The  objective  of  this  investigation,  therefore,  consists  of  deter- 
mining a system  of  partial  differential  equations  of  motion  for  thin 
elastic  bars  (or  foundations)  that  are  rigidly  clamped  on  an  entire 
lateral  surface.  Such  bars  will  hereinafter  be  called  "constrained 
bars."  These  equations  may  be  interpreted  as  describing  the  motion  of 
(i)  the  type  of  elastic  foundation  described  in  the  opening  paragraph 
of  this  section  or  (ii)  an  isotropic  crystal  bar  mounted  in  an  ultra- 
sonic device.  A variational  principle  is  used  to  obtain  the  field 
equations  for  the  system.  The  normal  and  shear  strain  components  of 
the  strain  energy  of  the  system  are  adjusted  by  means  of  the  introduc- 
tion of  correction  factors  of  the  type  that  are  used  in  the  Timoshenko 
and  Mindlin  theories  for  elastic  beams  and  plates,  respectively.  First 
order  and  second  order  theories  are  derived,  and  the  various  correction 
factors  are  evaluated  by  requiring  the  lowest  depth-shear  and  depth- 
stretch  frequencies,  as  derived  from  the  approximate  theories,  to  be 
identical  to  those  obtained  from  the  full  theory  of  generalized  plane 
stress  (see  Section  3).  Finally,  the  frequency  equation  for  a simply 
supported  constrained  bar  of  length  z is  derived,  and  plots  of  the 
variations  of  the  first  eight  frequencies  of  the  transverse  and  longi- 
tudinal modes  versus  the  length  to  depth  ratio  are  presented  and  dis- 
cussed. 
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2.  GENERALIZED  PLANE  STRESS 


Consider  the  bar  of  thickness  b in  the  x^  - direction  (b  <<  a , 
where  i denotes  a characteristic  length  of  the  bar)  that  is  depicted 
in  Figure  1.  The  bar  is  assumed  to  be  clamped  to  a rigid  surface  on 
the  face  = +h/2.  This  implies  that  the  longitudinal  and  transverse 
displacements  u,  (x-j , t)  and  u2  (x^ , x^,  t),  respectively,  are 
constrained  so  as  to  vanish  at  the  surface  x^  = h/2,  i.e., 

ua  (x-j , h/2,  t)  = 0,  a = 1 , 2.  (2) 

Furthermore,  the  lower  face  of  the  bar  is  free  of  applied  surface 
traction  and  its  longitudinal  and  transverse  displacements  are  com- 
pletely unconstrained.  Thus,  the  boundary  conditions  on  this  face  are 
those  of  vanishing  shear  stress  and  normal  stress,  namely. 


°12  = °22  = 0 on  x2  = ~h/2-  (3) 

The  equations  of  motion,  according  to  the  theory  of  generalized  plane- 
stress  [6],  for  the  in-plane  deformations  of  a thin,  elastic  plate  are 


a.=pu,  a,  B = 1 , 2 , 
aB»3  a 


(4) 


where  p denotes  the  density  of  the  material,  a . . = 3a  „/3xn,  and 

ap,p  an  d 

summation  over  repeated  indices  is  implied.  The  constitutive 
equation  is 


a = 2ye  + ye  6 0 , 
a&  aB  aa  aB 


where 


y = 2Xp/(A  + 2p) , 


(5) 


(6) 


^I.S.  Sokolnikoff.  Mathematical  Theory  of  Elasticity.  New  York: 
McGraw-Hill  Publishing  Co  (1956) . 
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Figure  1.  Coordinate  system  for  the  bar. 
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with  m and  x denoting  the  Lame  constants.  The  strains 
defined  by 


are 


e 


a6 


e«6  ' 1/2  + Va1' 


(7) 


It  is  easy  to  verify  that  the  above  field  equations  can  be 
derived  directly  from  a generalized  form  of  Hamilton's  principle  due 
to  Yu  [7],  (see  also  [8]),  linearized  and  specialized  here  for  the 
case  of  generalized  plane  stress.  This  principle  is  expressed  as 


6 / ' Ldt  = 0, 


where 

L = T-V+W 

is  the  generalized  Lagrangian  with 

T = J 1/2  pu‘aua  dA,  dA  = dx-jdx2, 


(°a»  Eae'  °aBeaB+  U)  dA> 


vj 

J A 

I p u ds  + p„(u  - u ) ds. 
C„  ra  a c otv  a a 


(8) 


(9) 

(10) 
(ID 

(12) 


Here  the  e „ are  generalized  strains  that  will  be  described  later, 

otg 

U,  defined  by 


Y.-Y.  Yu.  Generalized  Hamilton's  principle  and  variational  equation  of 
motion  in  nonlinear  elasticity  theory,  with  application  to  plate 
theory.  J.  Acoust.  Soc.  Amer.,  36.*  pp.  111-120  (1964). 

O 

G.  L.  Anderson.  The  stability  and  non-linear  vibrations  of  an 
orthotropic  hinged-hinged  beam.  Watervl iet  Arsenal  Technical  Report 
WVT-TR-75036  (1975). 
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(13) 


is  the  strain  energy  density,  which  is  assumed  to  exist  and  is  a func- 
tion of  e and  p = 0 n . where  n„  are  the  components  of  the  unit 

Otp*  Ot  Ot  p p p 

normal  vector  to  the  boundary  of  the  area  A.  Moreover,  t0  and  t^  are 
two  instants  of  time,  u = 8u  /at,  and  a barred  quantity  is  a prescribed 
function.  denotes  that  portion  of  the  boundary  curve  C on  which 
the  tractions  p are  prescribed,  Cu  that  portion  of  C on  which  the 
displacements  u^  are  prescribed,  and  s denotes  arc  length  along  C. 
Because  the  variations  eu^,  Se^,  are  arbitrary  throughout  the 
area  A of  the  body,  6u^  arbitrary  on  C^,  and  fip^  arbitrary  on  Cu,  their 
coefficients  in  the  varied  form  of  equation  (8)  must  vanish.  There- 
fore, the  present  form  of  Hamilton's  principle  will  generate  the  stress 
equations  of  motion,  the  stress-strain  relations,  the  strain-displace- 
ment relations,  the  stress  boundary  conditions,  and  the  displacement 
boundary  conditions.  As  will  become  apparent  in  subsequent  sections, 
this  process  will  be  quite  convenient  for  the  derivation  of  the  field 
equations  for  constrained  bars. 

3.  DEPTH  MODES  IN  AN  INFINITE  PLATE  STRIP 

The  equation  of  motion  in  equation  (4)  subject  to  the  boundary 
conditions  in  equations  (2)  and  (3)  admit  some  rather  elementary  solu- 
tions, herein  called  depth-shear  and  depth-stretch  modes,  in  the  case 
of  a plate  that  extends  to  positive  and  negative  infinity  in  the  x-j- 
direction.  The  lowest  frequencies  of  these  modes  will  be  of  consider- 
able importance  later  in  Section  5. 
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The  depth  modes  shall  be  characterized  by  a displacement  field 


that  is  independent  of  the  axial  coordinate  x,  i.e., 

I 9 

u = u (x?,  t).  (14) 

a a c. 

Therefore,  in  view  of  equations  (14),  (5),  and  (7),  the  stress 
equations  of  motion  in  equation  (4)  become  simply 

M^i  22  - p^i  | ^ ^ I (15) 

and 

C1 u2 , 22  = pU2’  |x2l<  h/2,  (16) 

where 

c-j  = 4u(u+x)/(2u+x),  (17) 

for  the  depth-shear  and  depth-stretch  modes,  respectively.  Moreover, 
by  virtue  of  equations  (2)  and  (3),  the  associated  boundary  conditions 
are  found  to  be 

u , = 0 on  x,  = -h/2,  u = 0 on  x,  = h/2  (18) 

a, 2 2a  2 

Assuming  solutions  in  the  form 

u (x0,  t)  = v (x,)  coswt,  (19) 

a 2 a c. 

where  u designates  the  natural  circular  frequency  of  vibration  of  the 
infinite,  constrained  plate  strip,  one  obtains  from  equations  (15), 
(16),  and  (18),  the  basic  boundary  value  problem 

v (x9)  + e2  v (x9)  = 0 (no  sum  on  a),  (20) 

0 t $ C.L.  L.  01  CL  L. 

V (-h/2)  = V (h/2)  = 0,  (21) 

Ot  9 Ot 

where 
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2 2 

^ = pu>  /p, 


(22) 


O _ 

^2  =po>  /c-j 


The  solution  of  equation  (20)  is 


v (x0)  = A cos3  x0  + B sin3x0. 

a 2 a a2  a a2 


(23) 


Substitution  of  this  expression  into  the  boundary  conditions  in 
equation  (21 ) yields 


\ 


A cos  3 h/2  + B sin  3 h/2  = 0, 

a a a a 

A sin  3 h/2  + B cos  3 h/2  = 0. 

a a a a 


(24) 


/ 

The  system  of  homogenous  equations  (24)  will  have  a non-trivial  solution 
if  and  only  if  the  coefficient  matrix  vanishes.  This  requirement  leads 
to 


cos3  h = 0, 

a 

whence 


(25) 


3 _ = (2n-l)n/2h,  n=l,2,3,... 

an 

Therefore,  from  equations  (22)  and  (26),  one  finds 


(26) 


n 2h  ’ p 


for  the  depth-shear  and  depth-stretch  modes,  respectively.  Of  particu- 
lar importance  later  will  be  the  lowest  depth-shear  frequency. 
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and  the  lowest  depth-stretch  frequency 


Finally,  it  may  be  noted  that 

A = sin  g h/2,  B = -cosg  h/2 

a a a a 

represents  a possible  solution  of  equation  (24).  Consequently, 
equation  (23)  becomes,  by  virtue  of  equation  (26), 

van  (x2)  = sin  [(2n-1)  * ( h— 2x2 )/4h] . (2 

Naturally,  these  are  the  eigenfunctions  of  equations  (20)  and  (21). 


4.  A THEORY  FOR  CONSTRAINED  BARS 


To  develop  a theory  for  constrained  bars,  one  can  proceed  with 
the  variational  principle  in  equation  (8)  in  much  the  same  way  as 
Vlasov  and  Leont'ev  did  in  their  book  [1].  Specifically,  it  is 
assumed  here  that  the  displacements  u^  can  be  approximated  by 


^1  (Xi  t)  f i (x^)  u (x-|,t)  + f 2 ( x 2 ^ ^ 

> 

u2(x1,x2,t)  = f-|(x2)  w (x i , t ) + f2(x2)  ^(Xpt), 


(30) 


Z.  Vlasov  and  l).  N.  Leont'ev.  Beams,  Plates  and  Shells  on  Elastic 
Foundations.  NASA  TT  P-  367  (19661 
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where  the  linearly  independent  coordinate  functions  f (x0)  are  re- 

a 2 

quired  to  satisfy  the  geometric  boundary  condition 

fa  (h/2)  =0,  a = 1,2,  (31) 


which  assures  that  equation  (2)  will  be  satisfied.  One  possible  set  of 
coordinate  functions,  whose  elements  satisfy  equation  (31),  is 

f 1 ( ^ 2 ) - 2h  * ^2^x2^  = 2 (^_^x2^^"*"^x2^ * (32) 

2h 

The  function  f-j  {x^)  has  previously  been  employed  in  references  [1]  to 
[4].  Of  course,  other  sets  of  functions  could  be  selected  and,  perhaps, 
would  be  even  more  desirable.  Further  comments  on  this  point  appear  in 
Section  9.  It  may  be  noted  that  the  functions  f (Xg)  given  in  equation 
(32)  are  orthogonal  over  the  interval  -h/2<X2<h/2,  i.e., 
f h/2 

J fi (x?)f?(x?)dx?  = 0 (33) 

-h/2  \ L c c c 

and  are  normalized  such  that  f (-h/2 )=1 . The  functions  f have  been 

a a 

plotted  as  functions  of  x = in  Figure  2.  Clearly,  f,  (Xg)  has  a 

single  zero,  whereas  f 2 (x£ ) *^as  two  zeros* 


^V.  Z.  Vlasov  and  U.  N.  Leont'ev.  Beams,  Plates  and  Shells  on  Elastic 
Foundations.  NASA  TT  F-357  (1 966). 

2 

N.S.V.  Kameswara  Rao,  Y.C.  Das,  and  M.  Anandakrishnan.  Variational 
approach  to  beams  on  elastic  foundations.  Proc.  Amer.  Soc.  Civ.  Engrs., 

J.  Engng.  Mech.  Div^97_,  pp.  271-294  (1971  ). 

3 

E.  H.  Dowell.  Dynamic  analysis  of  an  elastic  plate  on  a thin,  elastic 
foundation.  J.  Sound  Vib.,  35,  PP-  343-360  (1974). 

4 

N.S.V.  Kameswara  Rao,  Y.C.  Das,  and  M.  Anandakrishnan.  Dynamic  response 
of  beams  on  generalized  elastic  foundations.  Int.  J.  Solids  Struc- 
tures, 11_,  pp.  255-273  (1975). 
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Figure  2.  Variation  of  the  coordinate  functions  with  x. 


Substitution  of  equation  (30)  into  equation  (10)  leads  to  the 


following  expression  for  the  kinetic  energy: 

a 

T - \ ph  J [j  (02+w2)  + i (*2  +♦  2)]  dx, 


(34) 


By  equations  (7)  and  (30),  one  has 


ell  " fiu’i+V’l 


*22  fl  ,2  w + f2,2* 


'12  "*  2 ^1  2U  ^ 2,2^  ^l^’l  ^ * 1 ^ * 


so  that 


j aaBeaBdA  = J [u , -j N-j  + , -j + uQ^  + i|/Q2  + w , 1 Q3  + 4>,-|Q^  + 


+ wT-j  + <f>T2]  dx-j , 


(35) 


where  the  following  "generalized  stresses"  have  been  defined; 
h/2  _ h/2 


Ni  \ 1 ^o*  ~ \ f 1 1 * 

I _h/2  c L _h/2  C L 


h-  / 


h/2 


-h/2 


fl,2a12dx2’  Q2 


f h/2 

J ^2,2a~\ 2 dx2’ 


-h/2 


r h/2 

= J v 


1 2dx2’ 


-h/2 


/ 

-h/2 


h/2 


f 2^i 2^x2 * 


■,-/ 


h/2 


fl,2a22dx2’  X2 


-h/2 


I 

-h/2 


h/2 


f2.2  a22dx2‘ 
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Next,  suppose  that  one  defines  "generalized  strains" 


Y 1 1 , a]r  . . * , oi 2 2 such  that 

en  = Vll  + f 2 air  e22  = fl  ,2  y22  + f2,2a22’ 

2e12  = f 1 ,2  12  + f 2, 2 al  2 + flc12  + f2Cl  2* 

Then,  one  can  easily  show  that 

% 

J aaB6aBdA  = / ^UN1  +aHN2  + Y12Q1  + “l  2Q2  + C1 2Q3  + 

A 0 

+ ?12Q4  + y22T]  + a22T2]  dx] . 


(36) 


(37) 


Upon  writing  equation  (5)  as 


a „ = Z\ie  . + ye  6 „ 
ag  ag  aa  ag 


and  substituting  this  into  equation  (13),  it  follows  that 


2 2 2 
2U  - C-| e ii  + ^^2^11^22  ^ ^*1^22^  A^*^12* 

where  Ci  has  been  defined  in  equation  (17)  and 

c2  = 2yX/(X+2y) 


(38) 


(39) 


by  virtue  of  equation  (6).  In  order  to  correct  the  shear  deformation 
and  transverse  normal  strain,  correction  factors  <i  and  <2  are  intro- 
duced at  this  point.  Specifically,  e^  and  e22  shall  be  replaced  in 
equation  (38)  by  K^e^  and  K2e22,  respectively: 
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(40) 


2 2 2 2 2 
2U  = c-je^i  + 2K2C2e-|-|e22  +k2  cie22  + ^K1  lJe12' 

Inserting  equation  (36)  into  equation  (40)  and  integrating  the  result 
over  the  domain  -h/2<  x2  <h/2,  one  obtains 


2 2 7 

2U  dx2  ~ (h/3)Ciifi  i + (h/5)Ciai  1 _|C2c2Yl  1Y22~31<2C2Y1  la22 

1 2 2 
+ 3x2c2^r220lll  ~lC2c2al  la22+^^^K2  ^1^22 

+ (2/h)K22CiY22a22  + (l9/3h)K2^cla22+^^Kl^MY12+ 

+ ( 1 9/3h ) <i  ^uai  2+  (h/3)<i  ^y^i  2+(h/5)iCi  ^y^i  2 + 

2 2 12 

+ (2/h)K^  \iy^ 2al 2"k1  2^1 2+  3^1  ^12^12 

- 3^*1  yal 2^1 2 ya12^12‘  \4ly 


Therefore,  in  view  of  equations  (35),  (37),  and  (41),  equation 
(11)  becomes 
ri' 

V - J £u , i Ni  , i N 2^” u Q i +^2+*// 5 1 , i Q^"^wT-|  +(f>T ^ “ 

o 

- >i i Ni  -«i i N2_y1 2^1  ‘ “12^2  ‘ 2^3  “ C1 2^4  _ Y22T1 

~ a22^2  (^^^]Y']]  (^/^®)^]a]]  “ 2k2C2Y  1 1 Y22 


7 1 1 2 2 
~ ^ k 1 1 ot 22^"l6*<' 2^ 2^22^  1 1 "*  1^22^^ /2h)iC2  ^1*^22 

+ (l/h)K22CiY22a22  + (19/6h)K22cla22  + ^ /h)K]  2|JYl  2 + 
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12  2 2 
“2*1  l-1  Y -j 2^1 2 ^12^12  ~ (7/6)k.j  iJ0ti2^]2 


1 2 


(42) 


'^l  pa12?12  " 2K1 


Lastly,  because  tractions  may  be  applied  on  the  face  X2  = -h/2 
of  the  plate  and  because  of  equation  (2)  and  the  fact  that  f (-h/2)=l, 
equation  (12)  becomes  simply 


under  the  assumption  that  either  vanishing  tractions  or  displacements 
are  prescribed  on  the  ends  x-|  = 0,  £ of  the  plate. 

The  Lagrangian  function  L defined  in  equation  (9)  is  now  com- 
pletely known  because  the  required  expressions  for  T,  V,  and  W are 
given  in  equations  (34),  (42),  and  (43).  Therefore,  performing  the 
variational  process  indicated  in  equation  (8),  one  obtains  the  gen- 
eralized stress  equations  of  motion: 


which  are  subjected  to  one  boundary  condition  from  each  of  the  follow- 
ing sets: 


(43) 


Nl,l  - Q1  + P1  = 3 ph" 


N2,l  - Q2  + P1  = 5 ph^ 


(44) 


^3,1  “ T1  + P2  3 ph  W’ 
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(45) 


(i ) either  = 0 or  u = 0, 

(ii)  either  Ng  = 0 or  4>  = 0, 

(iii)  either  Q3  = 0 or  w = 0, 

(iv)  either  = 0 or  <{>  = 0 

at  x.j  = 0 and  =i  . The  generalized  stress-strain  relations  are 

^1  - ( h/ 3 ) c i "y 1 1 - 2 k2C2^22  " ^ K2^2a22’ 

N2  — ( h/ 5 ) c 1 1 + g*  k2c2Y22  ” 2K2c2a22 

Q-j  - ( 1 / h ) k 1 uy  12"*"  ( 2/  h ) k-|  ^“i  2 " 2 K1  ^^i 2 + g"  Ki  2’ 

Q2  = (19/3h)x.|2ua.|2  + (2/h)K12pY12  - (7/6)k.j  2p^  2 " ^<1 2^^i  2 ’ 

Q3  = (h/3)K1 2p c-j 2 " \ k]2my12  " (7/6)ic^wo^2»  (46) 

Q4  = (h/5)<-| 2uC-| 2 + (ki2/6)pY-|2  " (k-|2/2)uoi.|2> 

11  .2  2 
T 1 - — 2^2^2^11  ^ 6^  K2C2al  1 ^ (l/n)<2  c^  Y22  Cl / h ) K2 

7 1 2 2 

T3  — g-  k2c2y1  1 - 2 K2C2al  1 ^ 0/h)<2  ^l"^22  (19/3h)c^<2  *^22* 

and  the  generalized  strain-displacement  relations  are 

^11  = U,1  al  1 -^’l  Yi2  - U)  a] 2 ~ ^ ’ 


r22 


<*22  = $• 


Ifi 


c12  W’l 


512  **1 


w. 


(47) 


5.  THE  FIRST  ORDER  THEORY 


To  obtain  a first  order  theory,  one  sets  if>=  <p  = 0,  so  that  a „ = 

exp 

£>12  = 0 and  equations  (44)  to  (47)  reduce  to  (a)  the  generalized 
equations  of  motion 


Ni,i  - «i  + h ■ 3 

Q3,l  - T1  + P2  ■ J ph”’ 


0<X-|  <£ , 


(48) 


(b)  the  boundary  conditions 

(i ) either  = 0 or  u = 0, 

(ii)  either  = 0 or  w = 0, 

at  x-j  =0  and  x^  = a,  (c)  the  general 


(49) 

ized  stress-strain  relations 


N]  = (h/3)clYll 


2 k2C2y22! 


2 1 2 

Q"l  - (l/h)<i  py 1 2 ~2  K1  y?l 2, 

(50) 

Q3  = ( h/3 ) <i  ^C-|2"2k1  my12, 

1 2 
T-j  ~ " ~2  K2C2Y1 1 ciy22 


and  (d)  the  generalized  strain-displacement  relations 

Y11  ~ ^*1  Yi2  ~ ^ 2 - y22  - (51) 


Substitution  of  equation  (51)  into  equation  (50)  yields  expressions 
for  the  generalized  stresses  in  terms  of  the  generalized  displacements: 

N-j  — (h/3)  c -j u » i — 2 
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9 


(52) 


Tn  = - 77  K„C„U 


2 2 2 *1 


where  now  u(x^,t)  and  w(x.j,t)  may  be  interpreted  as  the  longitudinal 
and  transverse  deflections,  respectively,  of  the  face  x^  = -h/2  of  the 
bar.  The  displacement  form  of  the  equations  of  motion  can  next  be 
derived  simply  by  inserting  equation  (52)  into  equation  (48): 


To  evaluate  the  correction  factor  k-j  and  <2  in  equations  (53)  and 
(54),  one  possible  approach  consists  of  determining  the  lowest  natural 
frequencies  of  vibration  in  the  depth  modes  and  comparing  the  results 
with  equations  (27)  and  (28).  Hence,  if  one  assumes  that 

u = u(t),  w=w(t),  and  = 0,  then  equations  (53)  and  (54) 

become  simply 

-2  ..  p 

u + a)  u = 0,  w+u)W=0, 

where,  respectively, 

u>2  = 3K-|2y/ph2  and  w2  = 3K2^c-|/ph2. 

Comparison  of  these  expressions  with  equations  (27)  and  (28)  yields 


K-|2yw,ii-(3/h2)K22c-|W  + (3/2h)(<2c2-K^2y)u,^  + (3/h  )p£  =p'w.  (54) 


k j2  = < 22  = ir2/12. 


(55) 
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6.  THE  SECOND  ORDER  THEORY 


The  field  equations  for  the  second  order  theory  are  given  in 
equations  (44)  to  (47).  If  it  is  assumed  that,  for  an  infinitely  long 
bar,  u = u(t),  ...,<{>  =4> (t ) , then  the  deflection  forms  of  the  equations 
of  motion  are 


u + (3k.|  2y/ph2)  u + (6k i ^u/ph^ )i/j  = 0, 
(lOK^y/ph2)  u + $ + (95K.|2y/3ph2)^  = 0, 


(56) 


for  the  depth-shear  mode  and 

w + (SK^^c-j/ph2)  w + ( 3k^2Ci / ph 2 ) 4>  = 0, 

(57 

(5K22c1/Ph2)  w + * + (95K22Cl/3Ph2)  <f>  = 0, 

for  the  depth-stretch  mode.  If  solutions  of  equations  (56)  and  (57) 
are  sought  in  the  form 

u(t)  = A.|  cos  tot,  tp (t ) = B.|  cos  tot, 

w(t)  = A2  cos  tot,  <|>(t)  = B2  cos  tot, 


the  following  frequency  equations  are  obtained  in  the  usual  way: 

3(to/to0)4  - 104  (to/to0)2  + 105  = 0,  to02  “ic^y/ph2,  (58) 


for  the  depth-shear  modes  and 

3(to/toQ)4  - 1 04 (to/toQ ) ^ + 240  = 0,  toQ2  = K2^c-|/ph2,  (59) 

for  the  depth -stretch  inodes. 
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The  two  solutions  of  the  bi-quadratic  equation  (58)  are 

2 2 

K-pu 


2 < 1 u 

oj,  = 1.0408  ’ 

1 pfi 


Up  = 33.6251  pj 
P*i 


whereas,  according  to  the  expressions  derived  in  Section  3,  the  exact 
results  are 


J1 


p 2 

2 _ tt  y 


4ph 


2 ’ 


2 _ 9ir  y 

? 2 
4pli 


Thus,  one  might  require  that  the  values  of  be  identical,  so  that 
k-,2  =tt2/4  (1.0408)  and 


o^  - 


_ 5.684 


5V5- 


which  is  not  a good  approximation  to  the  exact  expression 


0)r 


For  the  depth-stretch  modes,  equation  (59)  yields 


k 2 

2 2 C1 
o)/  = 2.4859  — 2 


2 

2 k2  C1 
“2  - 32.1087 


whereas  the  exact  values  are 
p 2 c, 

2 _ IT  I 

W1  4 .2’ 

ph 


2 9*Z  C1 


2 4 ph2 


V V 

Proceeding  as  before,  one  finds  k ^ = tt  /4(2.4859)  and 


0)2  - 


3.594  tlVlL 

2 V .2 

ph 
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in  contrast  to  the  exact  value 


Clearly,  the  approximations  for  in  both  the  depth-shear  and 
depth-stretch  modes  are  not  as  accurate  as  might  be  desired.  Indeed, 
it  is  evident  from  equations  (56)  and  (57)  that  the  u,  \p  deflections 
are  coupled  in  the  depth-shear  modes  and  the  w,  <p  deflections  are 
coupled  in  the  depth-stretch  modes.  The  fact  is  clearly  at  variance 
with  the  corresponding  situation  as  embodied  in  equations  (15)  and  (16) 
obtained  from  the  "exact"  theory  of  generalized  plane  stress.  The 
source  of  this  difficulty  lies  in  the  choice  of  coordinate  functions 
f^x^)  made  in  equation  (32)  which  do  not  satisfy  the  orthogonality 
condition 

r h/2 

J fl ,2^x2^f2,2^x2^  dx2  = 0 (60) 

-h/2 


which  would  eliminate  all  the  undesirable  coupling  terms  in  equations 
(56)  and  (57).  Additional  comments  on  this  point  are  given  in  Section  9. 


7.  FREE  VIBRATIONS  OF  A CONSTRAINED  BAR 

2 2 

Suppose  that  p^  = p^  = 0 and  that  k(  = k , where,  according  to 
2 2 

equation  (55),  k = it  / 1 2 . Furthermore,  suppose  that  the  changes  of 
variables 

2 2 2 

X-j  = lx,  o<x<l , t = ct,  c = pfc  /<  ]1 

are  introduced  into  equations  (53)  and  (54).  Consequently,  the 
dimensionless  form  of  the  equations  of  free  motion  of  a constrained  bar 
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are,  according  to  the  first  order  theory, 

2 

3u"  - 3s  u + yw  = u, 
w"  - 3a?^w  - yiT  = w, 

where  u'  = 9u/9x,  u = 9u/9x,  etc.  and 

3 = cx/k2,  a = C-|/y,  5 = £/h,  y = (3i,/2h)  (<y-C2)/icy  (62) 

Suppose,  in  addition,  that  the  ends  = 0,  a of  the  beam  are  simply 
supported,  i.e., 

u = Q3  = 0 at  x-|  = 0,  A. 

In  view  of  equation  (52),  these  boundary  conditions  become  effectively 
u = w'  = 0 at  x = 0,  1 . (63) 

It  is  easily  verified  that  the  expressions 

u(x,x)  = A sin  n?rx  cos  tux,  w(x,x)  = B cos  nirx  cos  tot,  (64) 

where  n = 1,2,3,...,  A and  B are  arbitrary  constants,  and  to  is  the 
dimensionless  circular  frequency  parameter,  satisfy  both  the  boundary 
conditions  in  equation  (63).  Insertion  of  equation  (64)  into  equation 
(61)  yields  the  system  of  homogeneous  algebraic  equations 

[io2  - 3(nir)2  - 3£2]  A - niryB  = 0, 

-niryA  + [to2  - (n?r)2  - 3cx£2]  B = 0, 
which  has  non-trivial  solutions  if  and  only  if 
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-rm 


= 0. 


? 7 7 

w -g(n7r)  - 35 

2 2 2 

-n7ry  w - (mr)  - 3a5 

Expansion  of  this  determinant  yields  the  (dimensionless )frequency 
equation 

w4  - [(1+3)  (n7r)2  + 3(1  + a)  to2  + 3(nir)4  + 


+ [3(l+ag)52  - y2]  (mr)2  + 9a54  = 0, 
which  has  the  solutions 

(“j.n)2 ” 1 CO+eKn^)2  + 3(Ha){2]  + \ (-1 )J{[(l-e)(n„)2  - 

-3(1-c)52]2  + (2n„y)2  | 1/2,  j = 1,2,  (65) 

where  n and  n denote  the  (dimensionless)  natural  frequencies  of 
the  transverse  and  longitudinal  modes  of  vibration,  respectively.  By 
virtue  of  equations  (17),  (39),  and  (55),  the  quantities  defined  in 
equation  (62)  can  be  expressed  as 

a = 1/8  (1-v),  3 = 3/ 2tt2 ( 1 - v ) , y = (35/2)  [1  -W37ir(l  -v )]  , 
where  v denotes  Poisson's  ratio. 


For  large  values  of  the  1 ength-to-depth  ratio  5>  one  can  show  that 
2 

the  u>.  behave  asymptotically  as 
J 


03 


1 


2 

,n 


1- 


2v 

k(1-v) 


| +0(5’2) 


(66) 
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for  the  transverse  modes  and 


f 


^Z,n  = 3^2  + (niT)2<  3 + 4(1 -a) 


for  the  longitudinal  modes.  As  £+«>,  equations  (66)  and  (67)  become 


which  can  be  shown  to  be  the  dimensionless  forms  of  the  fundamental 
frequencies  of  the  depth-stretch  and  depth-shear  modes,  respectively, 
in  an  infinite,  constrained  plate  strip  (see  Section  3). 

In  Figures  3 and  4,  the  variations  of  u>.  have  been  plotted 

J » n 

against  £ for  v = 0.3  and  n = 1(1)8.  These  curves  show,  for  sufficiently 
large  values  of  £,  the  linear  rate  of  growth  in  frequency  as  indicated 
by  equation  (68).  For  j = 1 and  2,  the  frequency  of  the  lowest  mode, 
n = 1,  is  virtually  a straight  line  in  the  Sw-plane.  As  n increases, 
in  the  case  of  j=l,  the  frequency  curves  rise  steeply  with  increasing 
e,  and  then  abruptly  change  slope,  growing  thereafter  at  a slower  rate-. 

In  the  second  mode  (j=2),  the  higher  frequencies  increase  with  £ very 
slowly  until  an  abrupt  increase  in  slope  is  observed.  The  point  of  this 
transition  in  slope  occurs  for  larger  and  larger  values  of  5 as  the  mode 
number  n is  increased. 

8.  BEAMS  ON  AN  ELASTIC  FOUNDATION 

The  classical  equations  of  extensional  and  flexural  motion  of 
slender  elastic  beams  of  length  and  rectangular  cross  section  with 
depth  H and  thickness  b are,  respectively. 


“l,n  ' 5 ^ 


and  u>9  H /3  , 
£ j n 


(68) 
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£ 


Figure  3.  Variation  of  the  transverse  natural  frequency  w versus  f, 
for  the  firs  I eight  modes. 
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Figure  4 Variation  of  the  longitudinal  natural  frequency  versus  c 
for  the  first  eight  modes.  ^ 
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(69) 


and 


a2  'x2=+H/2’  A bH 


I = bH3/! 2,  E denotes 


Young's  modulus,  pb  the  mass  density  of  the  beam,  lT  (x^,t)  the 
longitudinal  deflection  and  w (x^,t)  the  transverse  deflection  of  the 
center  line  of  the  beam.  These  deflections  are  related  to  the  dis- 
placements u,  (x-|,x2,t)  and  u^Cx^.x^jt)  according  to 


u2  (x]  ,x2,t)  = w (x]  ,t), 
for  o<x^<^,  -H/2<x2<H/2. 

Under  the  hypothesis  that  the  beam  is  perfectly  bonded  to  the 
elastic  foundation  at  the  interface,  the  longitudinal  and  transverse 
deflections  must  be  equal  there,  i.e., 

IT  (x,,  H/2,  t)  = u (x,  , -h/2,  t), 

a I a I 

which,  in  view  of  equations  (30),  (32),  and  (71),  yield 


for  the  first  order  foundation  theory.  In  addition,  at  the  interface 


u1  (x.|  ,x2 , t ) = u (x.|  ,t)  - x2w51  (x^t). 


(71) 


u (x^ , t ) - (H/2)  w, i ( x i , t ) = u ( x i , t ) 
w (x,  ,t)  = w (x, ,t) 


(72) 


one  also  has  = -p  , where,  by  virtue  of  equations  (53)  and  (54), 

<2  ■< 
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(73) 


ai2  = (Ph/3)iJ'+(cih/3)u»1l  " (K^M/h)u+(K/2)(Ky-C2)w,-| 

(+ ) p 2 

a22  = " (ph/3)w+(K  yh/3)w,-| i_ (k  C2/h)w+(K/2)(c2-Kn)u,^ . (74) 

Some  special  cases  are  of  interest.  As  a first  example,  suppose 
that  the  beam  is  loaded  such  that  extensional  deformations  predominate 
over  flexural  deformations.  It  might  then  be  assumed  that  w = w = 0, 
so  that  the  pertinent  equations  of  motion  for  the  beam  and  the  founda- 
tion are  equation  (69)  and 

a-|2^  = - (ph/3)u  + (c-|h/3)u,j.|  - (<2u/h)  u.  (75) 

Substitution  of  equation  (75)  into  equation  (69)  and  use  of  lT(x^,t)  = 
u(x19t)  yield 

(EA+bhc^/3)u,-j i - (<2by/h)u  = (p^A  +pbh/3)ii  + ba^^  (76) 

As  a second  example,  suppose  that  the  beam  is  loaded  in  such  a 
manner  that  transverse  deformations  predominate  over  extensional 
deformations.  This  suggests  that,  as  a first  approximation,  u = u = 0 
and  w = w.  In  particular,  equation  (74)  may  be  approximated  by 

(+)  ..  2 2 
a22  = “(ph/3)^  + (K  ph/3 ) w,n  - (k  C2/h)w.  (77) 

Insertion  of  equation  (77)  into  equation  (70),  under  the  hypothesis  that 

(±) 

= 0,  leads  to 

EIw,.|iii  - (bhK2y/3)w,-| i + (bK2C2/h)w  + 

+ (p^A  + pbh/3)w  = - bc^^  ^ • (78) 
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Equation  (78)  possesses  a form  that  is  quite  similar  to  the  one 

investigated  in  references  [3],  [5],  and  [9]  to  [11].  As  cited  in 

equation  (1),  the  term  -(<  C2/h)w  in  equation  (77)  is  analogous  to  the 

Winkler  hypothesis,  whereas  the  set  of  terms  (<  uh/3)w,^  - (<  c?/h)w 

is  sometimes  identified  as  either  the  Wieghardt  or  Pasternak  model  for 

an  elastic  foundation.  If  the  inertia  term  - (ph/3)w  is  included  with 

these  latter  terms,  one  may  speak  of  a Wieghardt-  or  Pasternak-type 

inertial  foundation.  Essentially,  the  new  feature  in  equation  (78) 

2 2 

consists  of  the  introduction  of  the  correction  factor  k = n / 1 2 , which 
resulted  upon  the  insistence  that  the  foundation  model  should  produce 
the  fundamental  depth-stretch  frequency  for  an  infinitely  long  bar. 

In  the  general  coupled  case,  one  finds,  upon  inserting  equation 
(72)  into  equations  (73)  and  (74), 

a -(ph/3)u  + (c.h/3)u,ii  - ( <^u/h)u  + (phH/6)w,.  - 

12  i ii  i 

- (c  hH/6)wf111  + (</2h)  [<u(H+h)  -hc^w^,  (79) 
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a22(+)=  - (ph/3) w + (K/12)[Ky(4h+3H)-3HC2]w,11  - 

- (K^c2/h)w  + (k/2)(c2  - (80) 

Consequently,  substitution  of  equations  (79)  and  (80)  into  equations 
(69)  and  (70)  yields 

( EH+c-j h/3 ) u,n  - (K2y/h)u  - (c1hH/6)w,in  + 

+ (K/2h)[>y(H+h)-hc2]  w,^  = (p^H  +ph/3)  u - 

- (phH/6)  + a12  , (81) 

and 

- (c1bhH/6)u’,111  +(<b/2h)  [<y(H+h)  - hc2]  T7, ^ + 

+ (El  + c i bhH2/l 2)  w,1in  - (Kb/12h)[Ky(H+2h)2-4Hhc2]  w,n  + 

+ (<2c2b/h)  w = - (pbhH/6)  u,.,  - (pbbH+pbh/3)  w + 

+ (phbH2/l 2)  w,n  - b a22(_)  + ^bH  . (82) 

Neglecting  the  rotatory  inertia  terms  w,^  and  w,-^  as  well  as  u,^  in 
equations  (81)  and  (82),  one  obtains 

(EH  + c -j h/3 )TT , -j -j  - (K2y/h)  u - (c^H/Sjw,^  + 

+ (K/2h)[Ky(H+h)-hc2]w,1  = (pbH+phl3)u  + a12(_),  (83) 

and 

- (c.|bhH/6)u, -j -j i + (icb/2h)[icy(H+h)  - hc^lu,^  + 
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+ (EI+c1bhH2/l2)w,1111  - (Kb/12h)[<y(H+2h)2  - 4Hhc2]  w5]  + 

+ (K2c2b/h)w  = - (p^bH+pbh/3)  w -bo  ^ ^ bHo-|2  ^ *.  (84) 


9.  CONCLUSIONS 


In  Section  6,  it  was  observed  that  the  second  order  theory  does 
not  provide  accurate  approximations  for  the  second  depth-shear  and 
depth-stretch  frequencies  in  an  infinite  plate  strip.  This  difficulty 
was  traced  to  the  fact  that  the  selected  coordinate  functions  f (x„) 
stated  in  equation  (32)  do  not  satisfy  the  orthogonality  condition 
presented  in  equation  (60). 

The  undesirable  coupled  terms  in  the  equations  of  motion  can  be 

eliminated  upon  selecting  the  coordinate  functions  f (x0)  to  be 

\ 


f](x2)  = (h  -2x2)/2h, 

■^2^2^  = (^-2x2)  (h+^x2)  (h+10x2)/h2. 


(85) 


where  these  functions  satisfy  the  orthogonality  conditions  (33)  and 
(60).  If  a second  order  theory  is  required,  the  theory  can  be  reformu- 
lated following  the  procedure  outlined  in  Section  4,  provided  that  one 
now  use  the  coordinate  functions  given  in  equation  (85). 

On  the  other  hand,  instead  of  using  these  polynomial  coordinate 
functions,  it  seems  still  more  desirable  to  employ  the  eigenfunctions 
given  in  equation  (29)  as  the  coordinate  functions.  Specifically, 
higher  order  theories  for  constrained  bars  could  be  based  upon  expan- 
sions of  the  form 
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r-  /n) 

ua(xl  ,x2»t)  - / ( x i > t ) <{i  (X2) , 

n=l 

where,  in  view  of  equation  (29),  one  might  select 

<f>(n)(x2)  = H)n  sin  [ (2n-l ) Tr(h-2x2 )/4h] , -h/2<x2<h/2, 

as  coordinate  functions.  These  functions  satisfy  the  orthogonality 
conditions 


h/2 
■h/2 


/ dx0  = 0 if  n f k 


and 


/ 


h/2 
-h/2 


dx2  = 0 if  n / k. 
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